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We have investigated the time-evolution of a free particle in interaction with a phonon thermal 
bath, using the tight-binding approach. A dissipative quantum walk can be defined and many impor- 
tant non-equilibrium decoherence properties can be investigated analytically. The non-equilibrium 
statistics of a pure initial state have been studied. Our theoretical results indicate that the evolv- 
ing wave-packet shows the suppression of Anderson's boundaries (ballistic peaks) by the presence 
of dissipation. Many important relaxation properties can be studied quantitatively, such as von 
Neumann's entropy and quantum purity. In addition, we have studied Wigner's function. The 
time-dependent behavior of the quantum entanglement between a free particle -in the lattice- and 
the phonon bath has been characterized analytically. This result strongly suggests the non-trivial 
time-dependence of the off-diagonal elements of the reduced density matrix of the system. We have 
established a connection between the quantum decoherence and the dissipative parameter arising 
from interaction with the phonon bath. The time-dependent behavior of quantum correlations has 
also been pointed out, showing continuous transition from quantum random walk to classical random 
walk, when dissipation increases. 
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I. INTRODUCTION 

The nee to understand the evolution of a wave-packet (for instance, one particle in a one dimensional lattice) has 
been motivated by many problems in solid state physics pQ, quantum information [2], quantum open systems [3H6] 
and quantum optics [SJ [7], among others. In solid state physics in particular, the spread of a wave-packet from a 
highly localized initial state in the tight-binding lattice approximation [TJ [8] has drawn the attention of many authors 
[9j [10] . Interestingly, the concept of Quantum Walk (QW), borrowed from classical statistics [TTHT4] . has the same 
properties as a tight-binding free particle [TU [15]. Two kinds of QW are considered in the literature: discrete-time 
quantum coined walk [7J [9] [TTJ [I3j [T6H23] and continuous-time QW [TO] [I4j [I5j I24H27] . In the former (proposed by 
Aharonov et al.pj]), a two-level state, the so-called "coin", rules the unitary discrete-time evolution of a particle 
moving in a lattice. On the other hand, the evolution of the particle in the continuous-time QW is determined by a 
tight-binding like Hamiltonian [15] [27]. It is not difficult to see, by simple comparison, that a mapping between the 
tight-binding Hamiltonian and the QW model can be established (see Appendix A). 

There are some important differences between a Classical Random Walk (CRW) and a QW and these differences are 
well known in the literature [6j [9j [lOj [28] . The most important of these differences is the fact that for a tight-binding 
free particle the deviation <j(t) of the wave-packet becomes linear in time <j(t) ~ t (QW), in contrast to the CRW 
result that becomes similar to a(t) ~ \[t. This issue can be resolved straightforwardly, by analyzing, in the Heisenberg 
picture, the time evolution of the second moment of the position of the free particle [26] . QWs have been studied for 
possible applications in quantum information and quantum computation algorithms [2j [13] . Discrete-time evolution 
in random media has recently been used to study entanglement and quantum correlations in order to understand the 
role of noise and the mechanism of decoherence between the internal and spatial degrees of freedom in a QW [T7] . 

The study of a QW subjected to different sources of decoherence is an active topic that has been considered by 
several authors, in particular, due to their interest in understanding Laser cooling experiments [29 , modeling Blinking 
Statistics [30] . and carrying out quantum simulations [18]. Experimental study of quantum decoherence in discrete- 
time QW using single photons in the space was performed in [19]. These authors considered pure dephasing as a 
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decoherence mechanism and they could explore the quantum to classical transition by means of tunable decoherence. 
In other theoretical studies [lOj [I6j [20lf22] , the authors analyzed discrete-time and continuous-time QW in a random 
environment, and they could also study the quantum-classical transition. In [3Tj the authors study QW with decoher- 
ence by analyzing a non-unitary evolution in QW. Experimental analysis of QW in a random environment led to the 
study of quantum optical devices [jj [9] [23H25] ; interestingly, these experiments also show the non-classical behavior 
in DQW. 

In classical statistics and via the Central Limit Theorem, the Gaussian distribution plays a fundamental role in all 
random walks with finite mean square displacement per step [6j[28]. In quantum statistics, however, this analysis is 
much more complex because quantum thermal average must be taken using the reduced density matrix. On the other 
hand, in the Markov approximation, the time-evolution for the reduced density matrix requires a much more complex 
infinitesimal generator, and in fact this generator has been studied extensively, for many years, and is nowadays called 
the Kossakoski-Lindblad generator [32] [33] . One of the most interesting facts that distinguishes quantum mechanics 
from classical mechanics is the coherent superposition of distinct physical states. Many of the non-intuitive aspects of 
the quantum theory of matter can be traced to the coherent superposition feature. Related to this issue, an interesting 
question arises: How does coherent superposition operate in the presence of dissipation? These subjects have been 
important issues of research since the pioneer works of Feynman and Vernon [34] [35] and, Caldeira and Leggett [36] , 
among others, see for example the references cited in: [5] l6l I37U39] . 

An important fact in modeling a realistic QW is the inclusion of the quantum thermal bath B from the very 
beginning in order to get a dissipative open system [4-6 . To emphasize this fact we call this continuous-time model 
a Dissipative Quantum Walk (DQW), and this may be a mechanism of decoherence in a QW. Nevertheless there are 
many other mechanism of decoherence, see for example [31 . To our knowledge this pioneer problem was well posed 
in van Kampen's paper [12], but many other related approaches have also been presented in the literature [T5| 140]. 
The propagation of photons in waveguide lattices have been studied in recent years [24j [25 , and they are possible 
scenarios where our present results can be applied. 

In the present paper we study an open system in the Markov approximation to continuous-time [26] [41], i.e., a 
quantum mechanical particle that moves along a lattice by hopping while interacting with a thermal phonon bath £>. 
We have chosen an interaction Hamiltonian with the bath in such a way that this interaction produces the hopping in 
the tight-binding particle. Therefore, we highlight some of the issues of interpretation of the coherent superposition by 
tackling a soluble hopping model. The asymptotic long-time regime of the quantum probability, Wigner's function, 
quantum entropy, quantum purity, etc., are characterized as a function of the dissipation. The long-time decoherent 
behavior is also explained in terms of the present dissipative hopping model. This model analytically reproduces the 
(non-equilibrium) continuous transition from DQW to CRW when the diffusion coefficient goes to infinity (i.e.: when 
the temperature of the bath T — ^ oo). 



A. A tight-binding open model 

We have considered a free particle model (tight-binding approximation, for example see [8]) constrained to a one- 
dimensional regular and infinity lattice (in one band side) in interaction with a thermal bath of phonon B. The total 
Hamiltonian for this problem can be written in the form [T2] : 

/ i t\ 2 

H r = [E l-n^^j + % (1) 

The first term corresponds to the tight-binding Hamiltonian Hs where a and a) are translational operators in the 
Wannier bases \s) (it is easy to write these operators in second quantization, for more details see appendix [A| and 1 
is the identity operator. The second term is the interaction Hamiltonian and corresponds to a linear coupling between 
phonon operators B\ = B\ = VkBk and system operators V\ = = HTa, here T > is the coupling parameter. 

The third term is the phonon Hamiltonian written in terms of boson operators ^fouJkB^Bk [26] [41] [42] . Here Eq is 

k 

the tight-binding energy of site and Q the associated next neighbor hopping energy. 

The Quantum Master Equation (QME) for our DQW model can be obtained by eliminating the variables of the 
quantum thermal bath and assuming for the initial condition of the total density matrix Pt(0) = p(0) (g) p e % , where 
is the equilibrium density matrix of the quantum bath B. Then, in the Markov approximation, using [a, a"*"] = 
and a^a = 1, the evolution equation for the reduced density matrix is [26| 141] : 

p - % = ir t^//' d + D ( apat + a]pa ~ 2p ) ' ( 2 ) 
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with a trivial effective Hamiltonian: H e ff = Hs — hcu c l. The diffusion constant is given in terms of the quantum 
thermal bath temperature and the coupling constant in the form: 

D oc T 2 k B T/h, 

the additive energy hu c is related to the Caldeira and Legett frequency cut-off in the Ohmic approximation [36]. For 
simplicity we can add an additive constant to the tight-binding Hamiltonian — Eq + uu c h + ft. This assumption does 
not change the general results and finally we can write: 

Heff = n (1 - ^) , (3) 

as was presented in previous references p~2j |26l HT| . It can be noted from Eq.Q that as D — » (T 0), the von 
Neumann equation is recovered (unitary evolution). 



On £/ie second moment of the DQW 



From the QME Eq.Q, we can obtain the dynamics of any operator, in particular here we are interested in the 
evolution of the dispersion of the position operator q, which in the Wannier basis has the matrix elements (s is an 
eigenvalue of q): 

Hq|*'> = (4) 

note that q is defined as a dimensionless position operator with lattice parameter e = 1. Then the quantum thermal 
time-evolution of the first and second moments can be calculated straightforwardly. In the Heisenberg representation 
we get 

|q( t ) = ^ [q ,H eff] (5) 

|q 2 W = ^[q 2 ,iT e //]+2Z)l. (6) 

Where we have used that [q, H e ff] = ^ (a — a^) and that a(t) = a(0), a^(t) = a^(0), then we get for the time evolution 
of the position operator: 

qW = 9r(a-a f )* + q(o). 

It is simple to realize, just from a physical point of view, why there is no dissipative term in the equation of motion 
for q(t). Taking the thermal average in Eq. ^ we get ^ = ^ 5^s=^°oo ^ m [ps,s-i then introducing the 

explicit solution for the density matrix, Eq. ( [l3| ) , and using the properties of the Bessel function (with integer indices) 
^5=^00 Js+n (%) Js+m (x) = 5 n ,m- We get that ^ = indicating the conservation of reflection symmetry at 

the localized initial condition, otherwise if ^ (q(£)) 7^ 0, this would destroy the reflection symmetry principle. 

The quantum thermal statistical average -of any operator- in the Heisenberg picture can be written as (A(t)) = 
Tr [A(t)p(0)]. Then for the variance of the DQW we get: 

a(tf = <q(*) 2 > - (d(t)) 2 = ^ (y) 2 + 2Dt, (7) 

which is the expected dissipative result [6j [12j [26j [41]. From the Eqs. ([5| and ([6|, it is possible to see that von 
Neumann's term gives a contribution of the form oc t 2 for the time-evolution of the second moment, this is a well 
known quantum result. In fact for the null dissipation case, D = 0, we get that Anderson's boundaries (ballistic peak) 
movement is controlled by the linear law of the deviation of the wave-packet: 

a(t) = V /(q(t)2)-<q(t)) 2 = -l^, if D = 0. 

Therefore, we can associate the quantity Va = 4^ ^ with the velocity of Anderson's boundaries in a one-dimensional 
regular lattice. 



4 



In the appendix [C] we have calculated the second moment of the position operator by using the characteristic 
function, which is useful for calculating all quantum thermal moments. 

In the next section we will present, in detail, results concerning the probability profile of our quantum open model 
(D 7^ 0), i.e., a dissipative tight-binding free particle. We will study the time-evolution of the reduced density matrix 
[42j[43], characterize its decoherence, and solve, analytically, some correlation functions associated with the coherent 
superposition feature. 



II. TIME EVOLUTION OF THE DQW 
A. General properties 



In order to consider the evolution of our free particle in interaction with the quantum thermal bath £>, we have to 
solve the QME, Eq.(2|, for any time. This can be done in the Fourier representation. Let the Fourier "bra-kef 1 be 
defined in terms of the Wannier " bra-kef 1 in the form: 



1 oo 



-iks 



(s\. 



Then the QME adopts an explicit form: 



(ki\p\k 2 )= 



h 



(f fel -f fe2 )+2D(cos(fc 1 -fc 2 )-l) 



(ki\p\k 2 ) , 



(8) 



where = 0{1 — cos/c} . Note that the diagonal elements pk,k(t) = (fc|p(t) |fc) are constant in time, Pk,k(t) = 0, for 
example, for a localized initial condition p(0) = |sq)(so| (with sq = 0) we get: 



Pk,k(t) = pk,k(0) = vt, 

Z7T 



(9) 



even so, we can define a pseudo-momentum operator p = j^[q, H 3 ], where m represents the mass of the free particle 
in the model. We can calculate any moment of the pseudo-momentum operator, (p J ) = Tr '[p 3 p(t)] for j = 1,2, • • •. 
For the case j = 1, we obtain (p) = and for j = 2, we get (p 2 ) = ^(^) 2 m 2 . Then we can define the quantum 
thermal second moment of the velocity v in the following way: (v 2 ) = \(j-) 2 and so we re-obtain the velocity of 
Anderson's boundaries Va = ( see subsection "On the second moment of the DQW"). 
To solve Eq.® we define the function: 



F(k\, £2) = i-^ (cos(fei) — cos(/u2)) + 2D(cos(/ci — fe) — 1) , 
the general solution of the QME can be written as: 

(ki\p(t) \k 2 ) = p(0) kl k 2 exp(J r (/ci,/c 2 ) t) . 



(10) 



(11) 



In order to study the suppression of Anderson's boundaries by the presence of dissipation, it is convenient to go 
back to the Wannier representation (\s) = ^= dk e~ lks \k)). Adopting Eq.j9|, as initial condition for the density 
matrix, we get 



Hp(*)I*2> 



1 

2^ 



dk x / dk 2 e i{klSl - k2S ^ <fci|yo(*)|fc 2 > 



k 2 )(s 1 -S2)] p J 7 (k 1 ,k 2 )t 



We can solve Eq.(12) analytically if we consider Bessel's function properties: e lz 



E 



00 

n= — 00 



(12) 

i n J n (z)e ine ; 



Yln=-oo In( z ) e i where J n and /„ are Bessel functions of integer order. Using the following relations 
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J- n (x) = (-l) n J n (x), J n (-x) = (-l) n J n (x) and I- n (x) = I n (x), I n (-x) = (-l) n / n (x), where n is an integer 
[44j|45], we have obtained an analytical expression for (si\p(t)\s2) (for more details, see appendix |B|) : 



<S1|P(*)I*2> 



— »(S1-S 2 



n= — oo 



J s\-\-n 



nt 



' S2+n 



h J 



)I n (2Dt). 



(13) 



Using Xl^L-oo Jn( x ) — 1 an d X^^°oo -^ri(x) = e x [44 , we can check that Eq.(|13| fulfills the normalization condition 
Tr[p(t)] = ^2^L_ OQ (s\p(t)\s) = 1, VT>, and the probability of finding the particle in the site s in the lattice is: 



P s (t) = <s|p(t)|s> 



j, 



fit 



I„(2£>i) , DQW 



(14) 



It is straightforward to note that Eq.(13) contains all the information concerning the transition from DQW to CRW. 



In fact this transition is a genuine non-equilibrium one because any behavior characterizing the quantum to classical 
transition will be given in terms of the time-evolution of the density matrix, which of course is not a Gibbsian density 
matrix [3H5]. 

If D = (without dissipation, i.e., a closed system), we recover the QW and in this case the matrix elements 
(si\p(t)\s2) reduce to the form (using I n (0) = £ n? o, where n is an integer): 



Mp(*)I*2> d=0 = i {sl - S2) J Sl 



'fit 



fit 



, QW. 



(15) 



In the case D — >• oo, we can re-obtain exactly the classical probability [6j|28], where the off-diagonal elements of 
p(t) in Wannier basis are equal to zero. Alternatively, consider the limits Q — » and h — >• in such a way that: 

hm — 0, 



then from Eq.(14) and using J s+n (0) = J s+n5 o, where s + n is an integer it follows (for any finite time t) that 



- 2Dt E 

n= — oc 



[5 S> _„] 2 I n (2Dt) = e- 2I)t 



I s (2Dt) : 



which is just the probability of the CRW, i.e.,: 

P s (t) = e - 2Dt I s (2Dt) , CRW. 



(16) 



We note from Eq.dT6|), when t — » oo that we re-obtain the well known Gaussian asymptotic scaling for the CRW 



probability P s (t — >> oo) 



-1/2 



/4-kD 



(we have used the asymptotic limit of I n (x) 



= , for x oo [44 



Here we have considered it appropriate to define a new parameter such that tb = (rate of characteristic energy 



scales in the system) and t f = jj-t a dimensionless time, in order to plot the analytical expression (14). In Fig 1 w e 
show the probability at the site s P 3 (i) = (s\p(i)\s), for different values of the dissipative parameter (see, Eq.(^|). 
Similar plots for the probability profile in the presence of dissipation have been analyzed by Esposito et al. [T5] . 

Note from Fig {!] that for D = [rjj = 0), the system is not interacting with the phonon bath; in this case, we 
get a closed system, and therefore the evolution of the wave-packet is not diffusive but is ballistic, Eq.([l5]). So for 
t' = 31.8 we observe two maximum peaks in s p « ±29, which correspond to ballistic peaks (Anderson's boundaries), 
and far away from these peaks for |s| > s p the probability quickly goes to zero. In the case that |s| < s p , we easily 
see oscillatory behavior because of the quantum behavior of the system. If the dissipative term is small rjj « 1 
(td — 0.05) the oscillatory behavior starts to decrease and when the dissipation is of the order of the hopping energy 
(t\d = 1) or larger (ro = 5, 10, see Fig.l), the dissipation dominates in the system and the quantum character vanishes. 
In this case the wave-packet tends to a Gaussian form. This is the regimen for the CRW, Eq.(16). 



In the remainder of the paper we are going to study the probability profile of the DQW, the quantum purity, the 
Wigner function, and von Newmann's entropy as a function of the dissipative parameter vd- In the Wigner section 
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FIG. 1: Probability P s (t') for an initial condition p(0) — |so)(so| (where so = 0) as a function of position s, for m = 0, 0.1, 
0.5, 1, 5, 10 and t' — 31.8, where we have defined the new variables m — -^j^ and t' = j-t (dimensionless time). 



we are also going to introduce a criterion to describe the quantum-classical transition (see inset of figure 5). In a 
future work we will present analytical and numerical results such as the concurrence, negativity, etc., as a function of 
vd in order to study the decoherence and the entanglement in a bipartite system related to our DQW. 

In Figj2]we show the same facts as in Fig{T (for some values of rp) but in three dimensions (3D); i.e., we have 
included an extra axis for the time t' (this kind of graphic representation is usually called quantum carpet). We cut 
the axis of the probability for convenience (for example we don't show the probability for t' = 0, P 3 (0) = 1). In 
this quantum carpet we observe the transition from quantum regimen (Fig{l] (a)) to classical regimen (Figjl] (d)). 
Oscillatory behavior for small values of D (td « 1) is observed, and the oscillations in the probability start to 
disappear when is larger than one. 

On the other hand, the quantum purity (^q(^) = Tr[p(t) 2 ]), a quantity that provides information about whether 
the state is pure (Vq = 1) or mixed (Vq < 1), can be calculated analytically. In the present model quantum purity 
can be calculated using Bessel properties 



oo oo 



V Q (t) = Tv[p(t) 2 ]= ]T ]T ( Sl \p(t)\s 3 )(sMt)\si) 

S1 — — 00 ss = — oo 

1 f; J S1+ Y^V (2Dt) 
t f2j S3+n ^)j Sl+n (^y n '(2Dt). 

n t — — rvn ^ ' ^ ' 



i (si-s 3 ) e -2Df 



xi (s 3 - Sl ) e -2Dt- 

n' — — oo 

Using the relations Y^=-oo Jn+m(x)Jn+m>(x) = $m,m> and X^-oo J n( x ) = ^o(2x), we get 

V Q (t) = e- ADt I,(ADt). 
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The quantum purity takes the value one for D = (without dissipation) for all time, and for the case D ^ the 
quantum purity takes values smaller than one for t > (mixed state). For D ^ 0, the quantum purity takes on the 
asymptotic power law behavior Vq(t) ~ £ _1//2 , for t — » oo (using I n (x) ~ ^|==, for x — >• oo [44]), in agreement with 
previous results presented in the references [4T] . 

Another important measure that will indicate the presence of quantum behavior is the Wigner function. In the 
next subsection, we will study this function. 



B. The Wigner function 

The Wigner function was originally formulated as a quasiprobability for the position and momentum of the particle 
in quantum mechanics. For continuous variables X and K, representing position and momentum in the phase space, 
the Wigner function is given by [47j [48] : 

W(X,K,t) = ^ J dY(X + Y\p(t)\X — Y)e 2iKY , (17) 

where p(t) is the time-dependent density operator. The Wigner function could have negative values whereby it is 
considered a quasi joint probability of X and K, and the marginal distribution for X and K can be obtained in the 
usual form: 

J dKW(X,K,t) = (X\p(t)\X), (a) 

J dXW(X,K,t) = (K\p(t)\K), (b) (18) 
where (X\p(t)\X) and (K\p{t)\K) are the time-dependent probability densities for X and K. The normalization 



condition for p(t) can be checked from Eq.(17): 



J " dX J dKW(X,K,t) = Tr[p(t)} 



1. (19) 



In our DQW model the space is an infinite one-dimensional regular lattice. Thus, in this case, we propose the 
Wigner function as: 

oo 

W(s,k,t) = — J2 (s + s'\p(t)\s-s'y ks \ (20) 

s'= — oo 



where \s) is a Wannier basis. This definition fulfills the required conditions for the Wigner function, see Eqs. (18) 



and (19), using J2 S instead of J dX, where k G [— 7r,7r] (first Brillouin zone). 



Using Eq.([13j) in Eq.([20j), we can write the Wigner function as follows: 

p -2Dt co , , /Qf\ /Qf\ 

w(a, t) = ^— %2s e%ks Js + s 'Mt) Js - s 'A t m 2D£ > ' (21) 

n,s' — — oo ^ ^ ^ ^ 

then, after some algebra and using Bessel's properties [44U46] (Xl^L-oo eini Jn+m{x)J n {x) — J m (2x sin(7/2))e* /3m , 
where f3 = tt/2 — 7/2), we find: 



W(s,k,t) 



-2Dt 



2tt 



J2 J 2s+2n(2^sm^)l n (2Dt). 



h 



(22) 



Using Eqs.(18 a) and (18-b) in Eq.(22), we recover the probability densities for position s and momentum k 
respectively (see, eqs. (14) and (|9|). 

Interestingly, the Wigner function has information about the transition from DQW to CRW. Analyzing the case 
D = (without dissipation in the system) the Wigner function adopts the following form: 



W(s,k,t) QW 



1 

2^ 



2s 



Sit . k 

2— sm - 
ft 2 



D = 0. 



(23) 
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FIG. 2: (Color online) Representation in 3D of probability P s (t') for an initial condition p(0) = |so)(so| (where so = 0) as a 
function of position s and t\ for m = 0, 0.05, 1, 10 [(a)-(d)], where the blue regions indicate, approximately, value zero and 
red regions have high value for probability. The variables td and t' are the same as in Fig [l] 
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FIG. 3: (Color online) Normalized Wigner function W(s, k, t')/Wmax as a function of position s and momentum k, for m = 0, 
0.05, 1, 10 ([(a) - (d)]) and t' = 30, where W max is the maximum value of Wigner function for each value of m- We have 
considered m = and t' — jj-t (dimensionless) . The regions in blue indicate negative values, the ones in sky blue indicate, 
approximately, zero value and the red regions have the highest values of the Wigner function. 

In the pure diffusive regimen (Q = 0, i.e., the CRW), the Wigner function can be written as: 

„-2Dt 



W(s,k,t) C Rw 



2tt 



-I 3 (2Dt), ft = 0. 



Then from Eqs.(23) and (24), the Wigner function can be re- written as: 

00 

W(s,k,t) = 2n W(s + n,k,t) QW W(n,k,t) C Rw, 

n— — oo 

then, changing n — » — n and considering that I- n (x) = I n (x), we get 

00 

W{s,k,t) = 2n ^2 W(s-n,k,t) QW W(n,k,t) C Rw, 



(24) 



(25) 



(26) 



n= — 00 



this expression shows that the Wigner function of the DQW is given in terms of a non-trivial space convolution 
operation between the QW (D = 0) and the CRW (Q = 0). For the case D ^ this expression shows the nonlocality 
of the quantum mechanics of a free particle interacting with the quantum thermal B. 

In Fig{3) we show the Wigner function from Eq.([22|) or Eq.([26|. We have normalized the Wigner function with 
respect to its maximum value W max for each value of tb (i.e.: we use W(s, /c, t')/Wm ax , where t' = jj-t and tb = fn^)- 
Similar Wigner's quantum carpets have been analyzed for a QW on a ring without dissipation in [49]. In Figj3[a), 
we observe the quantum behavior for = 0. In this case the Wigner function has negative and positive values, and 
it shows remarkable oscillatory behavior (see Eq.(23)). In the case of small values of the dissipation, vd « 1, we 
observe (Figj3^b), with rjj = 0.05 << 1) that the oscillatory behavior starts to diminish and for values of > 1 (fig 
[3^c)-(d)) the dissipative term is dominant. In this case, the Wigner function takes only positive values, therefore it is 
a well-defined joint probability, and the system resembles a CRW when vd tends to infinity. In this case, the system 
is diffusive (see Eq.([24|)). 

In a complementary way, we now show the Wigner function (similar to Fig|3|, but we have fixed the parameter 
td = 0.05 and changed the dimensionless time t' . In Figj4[(a)-(d)], we observe the Wigner function for t' = 1, 10, 20, 
30. The time-evolution of the Wigner function also shows the transition from the DQW to the CRW. 
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FIG. 4: (Color online) Normalized Wigner function W(s, k, t')/Wmax as a function of position s and momentum k, for t' = 1, 
10, 20, 30 ([(a) — (d)]) and td — 0.05, in this case W m ax is the maximum value of Wigner function for each value of t' . 




FIG. 5: The Wigner function W = W(0,ir,lf), as a function of time t\ for r D = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 1, 2. The 
Wigner function W\ = W(0, 7r, 1.9), as a function of r^, see inset. For > 0.52 the Wigner function is non-negative, and well 
defined as a joint probability density in the phase space. 



In order to define a criterion to indicate when quantum correlations dominate over the classical correlations, or vice 
versa, we show in Fig{5]the Wigner function Wo = W(0,7r,t / ), as function of time t' . We note an oscillatory behavior 
of Wo for tb < 1- In the inset of Figj5j we show the Wigner function Wi = W(0,7r, 1.9) as a function of tb (we 
consider the value of t' where Wo takes their first minimum). We use this as a criterion to indicate when quantum 
correlations are more important than classical correlations, and from Figj5j we note that for tb ^ r< b-> with r c D = 0.52, 
the Wigner function is Wo < 0, therefore, the quantum correlations dominate in the system. For the case when 
tb > r c Dl the Wigner function is non-negative for all parameters in the system. In this case the classical correlations 
dominate over the quantum correlations. Therefore, in this situation, we could say that the Wigner function is a joint 
probability density in the phase space (because W(s, fc,£) > 0). This criterion is not unique, other criteria could be 
used. 

These results are consistent with the previous results obtained with the probability of finding the particle in the 
lattice site (see Figures [lp] ) . 
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C. The long-time limit of the reduced density matrix 



In this subsection, we have analyzed the long-time limit behavior of the reduced density matrix. The Eq.(12) can 
be written in the form: 



<*1|P(*)I*2> 



r dki ei ™ cos fclJ(fci)) 

J-7T 2tt 



(27) 



where 



X(fc 1 )=y" r ^e- i T cosfe2 cos[(fc 1 -A ;2 )( Sl - S2 )]e : 



2Dt(cos(fei-fe 2 )-l) 



(28) 



By symmetry the term proportional to sin [(fei — fc 2 ) (si — S2)] cancels out. Noting that in the long-time limit T{k\) 
can be calculated by using the stationary phase approximation [50] , we get (in the limit: Q/h » 2D and t — )> 00) 



X(fei) 



-2Dt(l- 



j fei)— i(fit— 7r/4) 



COS fciSi 



+e -2Dt(l+cos fcl )+^--/4) cos( ^ 



Si - 7T5 2 



>}■ 



Introducing this expression in Eq.(|27|) we can apply once again the stationary phase approximation (taking care that 
the saddle point (0,0) in the bi-dimensional integration does not contribute), then for the reduced density matrix in 
the asymptotic limit ((Q/h)t — >• 00) we get: 



2h 

(Sl\p(t)\s 2 ) -> ^^{cOS7r(5i +S 2 ) +COS7r(si -S 2 ) 



-4Dt 



[sin(2^£) cos(si + s 2 )^- cos(si - s 2 )^] 



-ze 



-4£>£ 



.0 



7T , 



[cos(2— t) sin(si + s 2 )- sin(si - ^2) ^D- 



(29) 



This result is in agreement with the asymptotic approximation of the Bessel function when it is replaced in the 

Eq.(13) (for x — >• 00, J n (x) ~ ^J~^cos(x — ^ — )). Thus we see that for long-time, terms proportional to e~ ADt will 

contribute to the quantum entropy production associated with the reduced density matrix. In fact, this contribution 
is proportional to the following matrix (in Wannier representation): 



1, 0, +!■ 



-1 

+1 



• A 


-iB 


A 


-iB 


A 


-iB 


A ■ 


■ iB 


A 


iB 


A 


iB 


A 


iB ■ 


■ A 


-iB 


A 


-iB 


A 


-iB 


A ■ 


■ iB 


A 


iB 


A 


iB 


A 


iB ■ 


■ A 


-iB 


A 


-iB 


A 


-iB 


A ■ 


■ iB 


A 


iB 


A 


iB 


A 


iB ■ 


■ A 


-iB 


A 


-iB 


A 


-iB 


A ■ 



(30) 



where A = ^(1 - sin(2^£//i)e- 4Dt ), B = ^ cos(2nt/h)e~ 4Dt . This means that asymptotically the probability 
profile p ss (t) goes to zero uniformly in the lattice, similar to ~ ^|^(1 — sin(2^t/^)e _4jD£ , and the off-diagonal elements 
form a time-dependent coherent binary structure of {A, ±iB} (in the Wannier representation) that also goes to zero 
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asymptotically as ~ 1/t. Due to the fact that Anderson's boundaries move at finite velocity Va = away from 

the initial condition (state \sq) (sq\), we expect that the DQW will be always inside a (time-dependent) finite domain 

of maximum size L ~ \J (q(t) 2 ) — (q(t)) 2 >> e (e is the lattice parameter, we take e = 1), which increases linearly 
in time. Then, we can calculate the asymptotic eigenvalues of (si\ p(t) ^2) by approximating p(t) to be a L x L 
finite-domain matrix. 

It is simple to calculate the non-null eigenvalues of a matrix of dimension L x L of the form ( |30| ). For D = 
(without dissipation) there is only one non-null eigenvalue: 

On the other hand, for the case D / 0, but fl/h » 2D, we obtain only two non-null eigenvalues for the reduced 
density matrix. These eigenvalues are as follows: 



±\/l + (L 2 -l)e- SDt -2Le- ADt sin(2^) + e - 8 - D * sin 2 (2^)}. 



(32) 



As expected, we observe from this expression that if we take D = 0, we recover the Eq.(31). These results allow us 
to calculate the long-time behavior of the von Neumann entropy. 

III. QUANTUM ENTANGLEMENT IN THE DQW 

A. Time evolution of von Neumann's Entropy 

In order to study the irreversibility behavior of a free particle in interaction with a quantum thermal bath 6, i.e., 
our DQW model, we have calculated von Neumann's entropy corresponding to non-equilibrium situations, assuming 
that the system is prepared in a highly localized state (i.e., we adopt the initial condition p(0) = \sq) (sq|, with so = 0). 

S(t) = -Ty[p(t)lnp(t)}. (33) 

In an open quantum system the reduced density matrix evolves in a Markovian approximation, following the QME 
(Eq.Q). Due to the dissipation, the reduced density matrix p(t) will not be diagonal at any time t > 0. The 
information of the quantum entanglement between the quantum thermal bath and our free particle -in the lattice- 



can be obtained from the reduced density matrix given in the Eq.(13). Due to the fact that the total system is a pure 
state, the von Neumann entropy for the reduced density matrix can be used to measure the entanglement. Then we 
should calculate von Neumann's entropy numerically using p(t) in the Wannier base. We have fixed the size of the 
chain to be L, with L « V^t where Va = ^7^? thus we have diagonalized the reduced density matrix. Then we can 
use the following expression for the quantum entropy: 

= In Ai, (34) 

i 

where A^ is an eigenvalue of the reduced density matrix (si\p(t)\s2) in the domain [— . We expect that for finite 
times (even in an infinite lattice) S(t) 7^ 00. We have carried out this calculation numerically and we have shown in 
Figl(3p(t) as a function of time t' = for different values of the dissipative parameter rp ~ 2D 



n/h' 

In Figj6ja), we observe the quantum entropy as function of tb and t' (3D visualization). As expected for tb = 
(D = 0), the quantum entropy is S(t') = for all t' > 0. In this case the quantum entanglement between the free 
particle in the lattice and the phonon bath is zero, which means we can write the wave function of the total system 
as a product of one state of the free particle in the lattice and one of the phonon bath (separable state). Another 
trivial result is the case that t' = 0, where the quantum entropy is zero for > (in this case p(t' = 0) = |so)(so|)- 
In the presence of dissipation the total system is no longer a separable state between the particle in the lattice 
and environment. The reduced density matrix p(t) is a mixed state (for D > 0), and the quantum entropy starts 
increasing in time for a fixed value of rp. We also note that for a fixed value of t\ the magnitude of the quantum 
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FIG. 6: (Color online) von Neumann's entropy as a function of t' and td (Fig (a)). This function measures the quantum 
entanglement between the particle in the lattice and the phono ns bath. We show a plot in 2D (Fig. (b)) of quantum entropy 
as a function of t' for different values of m — 0, 0.01, 0.05, 0.1, 0.5, 1, 2. 




10 20 30 40 



entropy increases as increases. The quantum entropy gives information about the transition from the DQW to the 
CRW, which is consistent with the results obtained above concerning the probability profile and the Wigner function. 
Finally, in Figj6^b), we show results of the quantum entropy as a function of t' and different values of td- 

Similar analyses to measure quantum correlations have been carried out by considering a free particle (in a lattice) 
with an additional internal degree of freedom: the quantum "coin". Thus, the correlation between the coin and the 
spatial degree of freedom have been studied in detail, and these models share some similarities with our results despite 
the fact that by considering the quantum coin as the thermal bath, the latter has a finite Hilbert space [T7JET]. 

Interestingly, by introducing the stationary phase approximation in Eq.(12), it is possible to obtain asymptotic 
behavior for («S2 1 p(t) then we can study the time-dependence of S(t) in the long-time regime analytically. 



1. The long-time limit of von Neumann's entropy 



Having the result of Eq.(29) we can re- write von Neumann's entropy in the long-time regime in the form: 

oo oo 



Si = — OO S2 = — 00 

L 



(35) 



here Xj are the eigenvalues associated with the asymptotic long-time regime of the reduced density matrix (si \ p(t) l^)- 
Using Eq.(32) in Eq.(35) and considering L » 1 (for simplicity, we re- normalized the eigenvalues of pit))] we get 
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(when § >> 2D and ft >> 1): 



-4Dt 



-4Dt 



2 

-ADt 



)ln(: 
1 - 



2 

-4L>£ 



-)ln( 



(36) 



In addition, another asymptotic approximation to S(t) can be made if we consider D ~ (with D£ << 1 ), so 



we can approximate e — »• 1 — 4Z}£ and replacing this expression in Eq.(36) we get a simpler expression for the 
quantum entropy in the form: 

S(t) ~ -2Dt\n(2Dt), Dt « 1. 

Note that in the limit D — >• (even when t ^> 1) the quantum entropy S(t) — » 0. We have checked this analytical 
result with our numerical calculation using Eq.(34) and the agreement is excellent, see Figlm 



IV. CONCLUSIONS 



A free particle -in an infinite regular lattice- in interaction with a thermal phonon bath has been studied by tracing 
out the degree of freedom of the bath. We then worked out the quantum master equation for a dissipative tight-binding 
model. 

We have solved the master equation analytically by using Bessel functions, and obtained the reduced density matrix 
p(t) as a function of the scale energies of the system (Q/H, D). We have also studied the transition from Dissipative 
Quantum Walk to Classical Random Walk in terms of parameters D and Q/h (or tb = fjjfi)' I n the case when 
2D « Q/H (td « 1) the quantum behavior is more important than the dissipation in the system. In the opposite 
case we re-obtain the Classical Random Walk 2D » Q/h (td » 1) because in this case the dissipation creates 
decoherence in the system. We have studied the Wigner function to analyze the pseudo-probability densities in the 
phase space. This function is very useful as an indicator of this quantum-classical transition. 

As an alternative approach to the study of the transition from Dissipative Quantum Walk to Classical Random 
Walk we have used tools from quantum information theory (as a function of dissipative parameter) to analyze the 
reduced density matrix. To describe this transition we have used von Neumann's entropy S(t) to measure the quantum 
entanglement between the free particle -in a lattice- and the phonon bath. We observed that for D = the quantum 
entropy is S(t') = for t' > (closed system), and when D increases we show that quantum decoherence starts 
to appear and therefore the quantum entropy increases in time with a law which is slower than that from classical 
statistics (S(t)cRw ~ m ^) [SI [28]. Asymptotically for D the quantum entropy turns out to be only a function of 
the dissipative parameter D. This fact also indicates the beginning of the transition from the Dissipative Quantum 
Walk to the Classical Random Walk. 

This analytical model allows us to study the effect of decoherence in the Dissipative Quantum Walk as a function 
of the two typical energies of the system. We can conclude that in the present model there are two characteristic 
time scales: the dissipative time td ~ 1/D and the hopping time th ~ h/Cl; the competition between these time 
scales controls the decoherence and correlation mechanism in the system. For example, the quantum purity Vq(t) is 
controlled by td, but in general the entropy and the interference phenomena appearing in the probability profile or 
in the Wigner phase- space pseudo-distribution are controlled by the competition between these time scales. 

The interesting problem of the propagation of photons in waveguide lattices are possible scenarios where our 
present results can be applied, also the analysis of the entanglement of a bipartite system can be studied in the 
present framework, works along these lines are in progress. In this way the present model gives insight into the effect 
of dissipation in more complex quantum systems, for instance, the analysis of quantum correlations between two 
particles -in a regular lattice- in interaction with a phonon bath. 
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acknowledges support received for this study from Universidad Nacional de Cuyo, Argentina, project SECTyP, and 
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Appendix A: On the second quantization and the one-particle tight-binding Hamiltonian 



A free Hamiltonian in the tight-binding approximation for spinless particles (fermion) can be written in second 
quantization in the form [52] : 

oo Q f 00 \ 

H S = E c l c * ~ 2 ( c l-i Cs + c Ui c s , ( A1 ) 

s = — 00 \s= — 00 / 

where c| and c s are creation and destruction operators in the site s of the lattice respectively (| • • • , 0, l s , 0, • • •) = cj|0), 
where 10) is the empty state). Then considering only one particle it is straightforward to compare Eq.(Al) with Hs 
in Eq.|l|, if we replace a and a) with a combination of c\, and <v, in the following way: 



a R = c l-i c s, = ^ c \+i c s, 



(A2) 



where c s are acting in the Fock-space. Therefore, we can also check that R and commute in the general case for 
many particles (where R^R = RR) = X^-oo c s c s — YlTs' -00 c l+i c I / -i c s c s / )? anc ^ f° r one particle in the lattice we 
get RE) = 1. 



Equation (A2) shows the expected mapping from Fock's space into the Winner basis. Then the connection between 



the tight-binding Hamiltonian and the QW model can be established. 



Appendix B: Reduced matrix density 



Here we show how to obtain Eq.(13). Replacing the following relations for the Bessel function: e 



iz cos 



J2n=-oc i n Jn(z)e in & ; e zcos6 = Xln=-oo ^i( z ) e in Eq.(12), where J n and I n are a Bessel functions of integer order 
[H|45], we find: 

(sMt)\s 2 ) = (^) e~ 2Dt £ J mi (^)j m2 (^V(2^) 

1 ,m2 ,n=— 00 ^ ' ^ ' 

jm 1 +m 2 I ^ e i(si+m2+n)fei / ^ —i(s2 — mi+n)k2 



Xl" 



Using the definition of Kronecker delta: 5 SjS ' = ^ e lk ^ s s ) in the previous expression, we finally obtain: 

00 / c 

( 5 i|pWk 2 )=^ (si - S2) e- 2 ^J Sl+ ^- 



Appendix C: Moments of the position operator - The characteristic function 

We have defined a characteristic function [6 for calculating moments of position operator q in the following way: 



g{0 = Tr\p(t)e* q ]= £ <*IP(*)l*>e*'> 



(CI) 



thus the quantum moments of q can be obtained using the following expression: 

1 <F 



<q(*) E 



MO 



Using Eq.(|13[) in Eq.( |Cl[ ), we can write the characteristic function in the form: 

ft 



g(£) = e -2r>t(i-oo8C) Jo ^ 2t 



— sin - 

h 2 



(C2) 



(C3) 
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here we have said that ^Z^ = _ 00 e m7 J n+m (x) J n (x) = J m (2x sin(7/2))e* /3m , where f3 = tt/2 — 7/2. From this char- 
acteristic function all the moments of the position operator can be calculated straightforwardly. In particular, we 
can re-obtain the variance of the DQW (see Eq.Q). Note that in the classical limit Q = we recover the expected 
characteristic function associated with the CRW j6J|28]. In general, Eq. (C3) shows that the non-equilibrium behavior 
of the characteristic function of the DQW is the product of the classical one and the quantum characteristic function 

Jo (2^ sin f). 
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